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Abstract Solutions to the BBGKY hierarchy of equations for molecular 
Brownian particle in ideal gas are considered, and exact relations are derived 
between probability distribution of path of the particle, its derivatives in re- 
spect to gas density and irreducible many-particle correlations of gas atoms 
with the path. It is shown that all the correlations always give equally im- 
portant contributions to evolution of the path distribution, and therefore the 
exact statistical mechanics theory does not reduce to classical kinetics even 
in the low-density limit. 



1 Introduction 

The idea of "molecular chaos" expressed by the Boltzmann's well-known 
"StoBzahlansatz" [l a is so much attractive that Bogolyubov, after he formu- 
lated [2] the exact hierarchy of evolution equations (now referred to as the 
Bogolyubov-Born-Green-Kirkwood-Yvon, or BBGKY, hierarchy [3]) for s- 
particle (s = 1,2,...) distribution functions, there and then truncated it 
at s = 2 to examine possibilities of substantiation of the Boltzmann equa- 
tion for dilute gases. In fact, however, until now the Boltzmann equation has 
no rigorous substantiation based on the BBGKy hierarchy, and role of the 
higher-order ( s > 2 ) distribution functions still is not properly understood!]] . 
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1 The frequently mentioned Lanford theorem [3] about gas of hard spheres under 
the Boltzmann-Grad limit concerns the so called "hard-sphere BBGKY hierarchy" 
(about it see also e.g. [3][MS]) which is not a true BBGKY hierarchy since represents 
interactions of the spheres by invented terms like the Boltzmann collision integrals 
(i.e. postulates what should be proved, if any). Besides, the Lanford's result spans 
too short evolution time intervals only. 
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In this communication I want to prove that no truncation of the BBGKY 
hierarchy can assert a qualitatively correct statistical description of gas evo- 
lution, even in case of arbitrary dilute gas. At that, in order to simplify the 
proof and at once essentially strengthen it, instead of the usual gas we will 
consider motion of a test particle in thermodynamically equilibrium ideal gas 
whose molecules interact with this particle only but not with each other (thus 
we concentrate on situation least favorable for inter-particle correlations). 

In Sec. 2 and Sec. 3 we formulate the BBGKY equations for this system 
and then rewrite them in terms of suitably defined "cumulant distribution 
functions" responsible for irreducible s-particle correlations between the test 
particle and s — 1 gas molecules ( s = 1 corresponds to the path probability 
distribution of the test particle). Next, show that these equations imply exact 
relations between any s-particle correlation (cumulant function) and deriva- 
tive of the previous (s — l)-particle one in respect to the gas density. In Sec. 4 
we will make sure that, consequently, the natural dimcnsionlcss measure of 
any of the correlations keeps non-zero even under the Boltzmann-Grad limit 
(hence, truncation of the BBGKY hierarchy always is incorrect). Finally, 
some other statistical properties of the correlations will be discussed. 



2 The BBGKY equations and cumulant distribution functions 

We want to consider thermal random motion of a test molecule (TM) in 
thermodynamically equilibrium gas, with specified position R(t) of the TM 
at some initial time moment t — : R(0) = Ro • 

Let P and M denote momentum and mass of TM, m , Vj and pj 
(j = 1,2,...) denote masses, coordinates and momenta of other molecules, 
^(r) is (short-range repulsive) potential of interaction between any of them 
and TM, and n is gas density (mean concentration of molecules). At ar- 
bitrary time t > , full statistical description of this system is presented 
by the chain of (k + l)-particle distribution functions (k — 0,1,2,...): 
Fo(t, R, P| Ro ; n ) which is normalized (to unit) density of probability distri- 
bution of TM's variables, and F k (t, R, rW , P, pW | R ; n) (where r-W = 
{ri... r,t } , p( fc ) = {pi... Pk } ) which are probability densities of finding TM 
at point R with momentum P and simultaneously finding out some k 
molecules at points Tj with momenta pj . A rigorous definition of such dis- 
tribution functions (DF) was done in 0. In respect to the coordinates rj they 
are not normalized, but instead (as in [2]) obey the conditions of decoupling 
of inter-particle correlations under spatial separation of particles (in other 
words, DF satisfy a cluster property with respect to spacial variables). Sub- 
ject to the symmetry of DF in respect to xj = {rj, pj} these conditions can 
be compactly written as follows: — > Fk-i G m (Pfc) at r& — ► oo , where 
G m (p) is the Maxwell momentum distribution of a particle with mass m . 

The enumerated DF satisfy a standard chain of equations [2]: 

¥ = [4fi, + 4|f (R_%l)Fw ' (1) 
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with k = 0, 1, ... and along with obvious initial conditions 

F k \ t=a =<5(R-Ro) exp(-.ff fe /T) = 
= 5(R - Rq) G M (P) nj=i S fo - R ) G ™(Pj) > 

where TJfc is Hamiltonian of subsystem " k molecules + TM " , [...,...] means 
the Poisson brackets, J, ... — J J ... dvk dpt , #'(r) = V#(r) , and _E(r) = 
exp [— <P(r)/T] . Notice that TM can be considered as a molecule of non- 
uniformly distributed impurity, and equations ((T|) are identical to the equa- 
tions of two-component gas |2j in the limit of infinitely rare impurity, when 
the main component is in spatially homogeneous and thermodynamically 
equilibrium state. 

Equations (JTJ) together with Q unambiguously determine evolution of Fq 
and eventually the probability distribution of TM's displacement R — Ro . 
These equations will become more clear if we make a linear change of DF 
-F/c by new functions 14 with the help of recurrent relations as follow: 

F„(t,R,P|Ro;n) = F (t,R,P| R ;n) , 

Pi(i,R,n,P,p 1 |Ro5n) = (3) 
= Vo(*,R,P|Ro;n)/(ri-R,pi) + V x {t, R, n, P, P i| Ro; n) , 

where /(r, p) = E(r) G m (p) , 

F 2 (t, R, rW, P, p« |R ; n) = V Q (t, R, P|R ; n) f((n,pi) /(p 2 , p 2 ) + 
+ Vi(t, R, n, P, pi|Ro; n) /(pa, Pa) + Vi(t, R, r 2 , P, p 2 |R ; n) f(p h Pl ) + 
+ y 2 (i,R,r( 2 ),P,p( 2 )|R ;n) , 

where pj = rj ~ R , and so on. 

Apparently, from viewpoint of the probability theory, Vk represent a kind 
of cumulants, or semi- invariants, or cumulant functions (CF). It is important 
to notice that if all these CF were zeros then all conditional DF of gas, 
Fk/Fo , would be independent on initial position Ro of TM and thus on its 
displacement R — Ro . This fact makes visible very interesting specificity of 
the CF Vk : they are not mere correlations between instant dynamic states 
of TM and k gas molecules but simultaneously their irreducible correlations 
with the total past TM's displacement. 



3 Evolution of many-particle correlations and their relation to 
density derivatives of the path probability distribution 

In terms of the CF the BBGKY hierarchy acquires a more complicated tridi- 
agonal structure (we omit uninteresting algebraic details): 

^ = [H h , V k ] + n A J k+ <P'(R - v k+1 )V k+1 + 
fe 

+ T]T P kj G m (pk)E'(rk-TL) 

3=1 



P _d_ 



Vk- 



(4) 
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Here E'(r) = VE(r) , and P k j symbolizes transposition of the pairs of 
arguments Xj and x k ■ On the other hand, initial conditions ([2|) and the 
above-mentioned clustering conditions [2 take very simple form: 

K,(0,R,P|R o ; n) =<5(R-Ro)G M (P) , 

y fc (0,R,r«,P,p( fc )|R ;n) =0 , (5) 
yfc(i,R,r( fc ),P,pW|R ;n) -> at v 3 , -> oo 

( 1 < j < fc). Thus, as it should be with cumulants, CF 14 disappear under 
removal of already one of molecules. 

From equations (j4|) as combined with the boundary and initial conditions 
([5]) it is clear that passage to the limit in ([5]) realizes in an integrable way, so 
that integrals V k = J k+1 Vk+i are finite. Let us consider them. By applying 
the operation J k to equations one easy obtains equations 

^ = [H k , V k ] + n A jf £'(R _ rfe+1 ) y k+1 + 

+ ^J^'(B.-T k+1 )V k+1 + (6) 

k 

+ PkjG m (p k ) E'(r k -n) 

3=1 

(with fc = 0, 1, ... ). Because of ([5]) initial conditions to these equations are 
zero: V k {t = 0) = at any k . 

Now, in addition to Vk , let us consider derivatives of CF in respect to 
the gas density, V k ' = dV k /dn . It is easy to see that differentiation of 
in respect to n yields equations for the V k ' which exactly coincide with ([6j) 
after changing there V k by Vl . Besides, in view of ([5]), initial conditions 
to these equations again all are zero: V k (t = 0) = at any k > . These 
observations strictly imply exact equalities V k ' — V k , or 

^- y fe (*,R,r«,P,p( fe )|R ;n) = (7) 

= / ^+i(t ) R,r(*+ 1 ) ) P,p( fc + 1 )|R ;n) . 

Jk+l 

This is main formal result of the present paper. Notice that it evidently con- 
firms the assumed finiteness of integrals V k = J k+1 V k +i . Together with the 
CF's definition ([3]), it forms a lemma on the way to similar general theorems 
of statistical kinetics of fluids. 

4 Discussion and resume 

The equalities (O contain the proof promised in Sec.l. Indeed, they show, 
firstly, that all the many-particle correlations between gas molecules and 



MT 



d_ 
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Vk- 



5 



total path, or displacement, of the test molecule (TM) really exist, i.e. differ 
from zero. Secondly, integral values of all the correlations, represented in 
the natural dimensionless form, have roughly one and the same order of 
magnitude. Indeed, multiplying equalities by n k and integrating them 
over TM's momentum and all gas variables, we have 

n k V k (t,A;n) = n k J...JJv k dP = n k ~ c k V Q (t,A;n) , 

where Vo(t,A;n) = J Vo(t, R, P|Ro; n) dP is just the probability density 
distribution of the TM's displacement, A — R — Ro , and c k some numeric 
coefficients obviously comparable with unit. Equivalently, unifying all CF 
into one generating function, we can write 

oo k k 

V (t, A-r^ + Y, 1 ^- V k (t, A; n) = V (t, A: (1 + u) n) . (8) 

fe=i 

Thus, distribution of total of the TM's random walk "is made of its correla- 
tions" with gas molecules like the walk itself is made of collisions with them, 
and hardly some of the correlations can be neglected if we aim at completely 
correct analysis of the BBGKY equations. 

We see also that characteristic volume occupied by the correlations has 
an order of the specific volume: ( | H ... J J Vk dP \/Vo ) x ^ k ~ n _1 . In the 
Boltzmann-Grad limit, n — > oo , ro — > , ixr^n = 1/A = const , it becomes 
vanishingly small as measured by the TM's mean free path A . But, never- 
theless, it remains on order of effective volume of the "collision cylinder", 
~ irr^X. This observation prompts that fc-particle correlations are concen- 
trated just at those particular subsets of fc-particle phase space which cor- 
respond to (real or virtual) collisions, and therefore their action holds under 
the limit. The same is said by the equality ([8]) which also holds out. This 
becomes obvious if we take into account that actually important parameter 
of the integrated CF under the Boltzmann-Grad limit must be A instead of 
n and rewrite © in the form 

OO k 

W (t, A;X)+J2^ W k (t, A; A) = W (t, A; A/(l + u)) , 

where Wk(t, A; A) = lim n h Vk(t, A; n) . Thus, nothing changes under the 
Boltzmann-Grad limit. 

It is necessary to underline that the correlations under our attention are 
qualitatively different from correlations which appear in standard approxi- 
mations of the BBGKY hierarchy and connect velocities of particles after 
collision (see e.g. [7])- In our notations, a pair correlation of such the kind 
would look nearly as Vi(t, R, r, P, p) = F 1 (t 7 R', r', P', p')—F 1 (t, R, r, P, p) , 
where the primed variables describe the pre-collision state. Clearly, because 
of the phase volume conservation during collision, integration of this expres- 
sion over p — v — R and the momenta results in zero. This observation shows 
that our correlations do live in the configurational space and connect coordi- 
nates and walks of particles (may be coexisting with statistical independence 
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of their velocities). The role and statistical meaning of such correlations were 
under investigation already in [8] (and in principle even earlier in [9]). By 
their nature, they are attributes of spatially inhomogeneous states and evo- 
lutions (evolution of Vo(t,A;n) gives an example). 

In view of all the aforesaid, we can suppose that the Boltzmann-Lorentz 
equation 3,5 and, moreover, the Boltzmann equation in itself and its gen- 
eralizations do not represent a (low-density) limit of the exact statistical 
mechanical theory. The classical kinetics is only its simplified probabilistic 
model (may be excellent in one respects and caricature in others). Of course, 
in the exact theory also molecular chaos does prevail. But here it is much 
more rich, even in dilute gas, and does not keep within naive probabilistic 
schemes. 



Appendix 

In the following, I append the Introduction and Discussion and resume sections 
from the first variant of this paper (titled "A truth of molecular chaos") rejected 
by the JSP (without any review or explanations) and by the CMP, on the grounds 
of that "it is in contrast with Lanford's theorem" and it "is not clear because 
assumptions, heuristic ideas and non rigorous steps are not clearly distinguished" . 

In order to conform to the latter remark, now I moved off all the heuristics 
which was localized in the Introduction and Discussion and resume sections (other 
sections remain exactly as before). Nevertheless, in my opinion, this heuristics is 
useful, and readers can see it below. What is for the Lanford's theore m, i t was 
sharply commented in the above footnote. For more detail comments see [10]) and 
also [8] where principles of correct collisional (with the help of the Boltzamnn 
integrals) description of spatially inhomogeneous evolution were formulated. Here, 
I can add that all the results of the present paper easy extend to hard-sphere 
interaction. 



Introduction 

One of creators of the modern probability theory A. Kolmogorov underscored [11] 
that in it " the concept of independence of experiments fills most important place " 
and " correspondingly one of most important objectives of philosophy of natural 
sciences " is " clearing-up and refinement of those prerequisites under which one can 
treat given phenomena as independent''^. Recall that in the probability theory some 
random phenomena or quantities A and B by definition are termed "independent" 
if th eir probability distributions are independent, that is P(A, B) — P(A) P(B) 
[llj . But in natural sciences the independence of phenomena A and B is thought 
as absence of cause- and- effect connections between them, that is absence of their 
mutual influence. Does independence in this natural sense mean independence in 
the sense of the probability theory? 

Certainly not from the viewpoints of common sense and philosophy. Merely 
because A and B which do not directly influence one on another nevertheless 
both can be parts of the same random event and thus turn out to be indirectly 
connected. 

From the scientific point of view, it is natural to bring the same question to the 
statistical mechanics. One of creators of modern theor y of dynamical systems and 
statistical mechanics N. Krylov thoroughly analyzed it [T2] and confirmed the neg- 
ative answer. He concluded that opinions that " phenomena which are "manifestly 



2 The italics in quotes means citations freely translated by me from Russian texts. 
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independent" should have independent probability distributions" , and the like, are 
nothing but " prejudices" |12j . 

Especially Krylov pointed [12] to the firmness of such prejudices^. Only it ex- 
plains why the molecular chaos hypothesis put forward by Boltzmaim many years 
ag o [1 ] until now dominates kinetics although never was somehow substantiated 
3,13 . And why N. Bogolyubov, when he obtained [2] an exact hierarchy of evo- 
lution equations for s -particle distribution functions, straight away truncated his 
equations at s — 2 taking in mind their reduction to the Boltzmann equation. 

Undoubtedly, molecules of rarefied gas are independent in the natural sense 
since almost surely have nothing common in the past. Nevertheless they can be es- 
sentially dependent in the sense of the probability theory. This is quite understand- 
able [8]. As colliding particles have no common history, there is no back reaction 
of the gas to past collisions of any of them. Therefore arbitrary long fluctuations 
in relative frequency of collisions are allowable^- These fluctuations just play the 
role of aforesaid random events producing indirect statistical interdependencies be- 
tween pairs (or groups) of particles capable of being participators of one and the 
same collision (or a cluster of successive collisions). 

As the consequence, P(A, B) ^ P(A) P(B) where P(A) is probability of find- 
ing a molecule at (phase) point A and P(A, B) is probability of finding simultane- 
ously two molecules at points A and B . At that, relaxation of one-particle distri- 
bution P(A) is determined by the two-particle correlation P(A, B) — P(A) P(B) . 
Relaxation of the latter just similarly always (regardless of the gas rarefaction) is 
determined by three-particle correlation. And so on up to infinity. Since during time 
interval t a molecule undergoes ~ t/r collisions (with r being characteristic free- 
flight time) , a correct description of gas evolution over this interval requires taking 
into account s-particle correlations with at least s < t/r . Hence, in practice, in 
contrast to the conventional opinion, the whole hierarchy of equations deduced by 
Bogolyubov [2] is necessary. 

In [8] and other worksQ approximate solutions to this hierarchy or, in other 
words, the Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) equations [3] were 
suggested for the problem about random wandering of a test molecule, and expla- 
nations were expounded why the Boltzmann's hypothesis is wrong. The aim of the 
present communication is to prove the statements of preceding paragraph without 
any approximations. At that we will strengthen the proof and besides simplify it 
due to replacing the usual gas by ideal gas whose molecules interact with the test 
molecule only but not with each other. 



Discussion and resume 

The result (0 contains the proof promised in Sec.l. Indeed, equalities ((7]l show, 
firstly, that all the many-particle correlations between gas molecules and past dis- 
placement of test molecule (TM) really exist, i.e. differ from zero. Secondly, all they 
have roughly one and the same order of magnitude. For instance, if comparing their 
integral values, due to (0 we can write, in natural dimensionless units, 

n j j jv k dP = n k V Q W (t,A;n) ~ c k V (t, A; n) , 

3 That "are so habitual that even persons who agreed with my argumentation 
usually automatically go back to them when facing with a new question" . 

4 "... relative frequencies of some phenomenon along a given phase trajectory, 
generally speaking, in no way are connected to probabilities" [12] , 

5 Kuzovlev Yu. E. : On statistics and 1/f noise of Brownian motion in 
Boltzmann-Grad gas and finite gas on torus. Part I. Infinite g as, arXiv: 
|cond-mat / 06095 1 5 ; Part II. Finite gas, arXiv: |cond-mat/ 06 12325 . See also 
Kuzovlev Yu. 1H. : Kinetic theory beyond conventional approximations and 1/f- 
noise, arXiv: |cond-mat /9903350 1 . 
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where V (t, A;ri) = / V (t, R, P|R ; n) dP is probability distribution of the TM's 

displacement Zi = R-R , V^ k) (t, A; n) = d k Vo{t, A; n)/dn k are its derivatives 
in respect to gas density n , and Ck some numeric coefficients. Hence, all the 
correlations are equally important, and none of them can be neglected if we aim at 
knowledge about true statistics of TM's random walk. 

For more details, let us suppose that (s + l)-particle correlation is so insignifi- 
cant that one can assign V s = in (j4]). At that, according to (0])-([5|, all higher-order 
correlations also will be rejected. Then, obviously, according to @, distribution 
Vo(t, R, P]Ro; n) and thus Vo(t, A;n) must depend on n definitely as an (s — 1)- 
order polynomial. But, from the other hand, distribution Vb what follows from 
the truncated chain of equations ([4]) certainly is absolutely non-polynomial func- 
tion of n . With taking into account that equalities (0 do express exact properties 
of solutions to Q-© we see that very deep contradiction is on hand. 

This contradiction clearly prompts us that truncation of the BBGKY hierarchy 
leads to qualitative losses in its solution. 

Some possible losses already were characterized in [8] (and principally even 
much earlier in [51ll4j) and in part filled up in [S] as well as preprints 4 (it may be 
useful also to see some of my recent worksQ). Therefore here we confine ourselves 
(continuing 5-th paragraph of Introduction) by remark that cutting of the (s + 
l)-particle correlation means cutting of s-th and higher statistical moments of 
fluctuations in relative frequency of TM's collisions with gas molecules (in other 
words, fluctuations in diffusivity of TM [9]). At s = 2 these fluctuations are 
completely ignored, and such truncated equations Q yield a closed equation for 
Vo(t, R, P|Ro; n) which is equivalent to the Boltzmann-Lorentz equation [3]. 

It is necessary to emphasize that above reasonings, as well as the exact relations 
(O, are indifferent to a degree of the gas rarefaction. Consequently, one can state 
that the Boltzmann-Lorentz equation (moreover, all the classical kinetics including 
the Boltzmann equation and its generalizations) does not represent a (low-density) 
limit of the exact statistical mechanical theory. The conventional kinetics is only 
(more or less adequate or caricature) probabilistic model of exact theory. Of course, 
in the latter also molecular chaos does prevail. But here it is much more rich, even 
if speaking about rarefied gas, and does not keep within naive probabilistic logics. 



6 Kuzovlev Yu. E. : Virial expansion of molecular Brownian motion versus tales 
of statistical independency, arXiv: 0802. 0288 ; Thermodynamic restrictions on 
statistics of molecular random walks, arXiv: 0803.0301 ; Molecular random walks 
in a fluid and an invariance group of the Bogolyubov generating functional equa- 
tion, arXiv: 0804.2023 ; Molecular random walks and invariance group of the Bo- 
golyubov equation (to appear) . In these works (see also references therein) another 
ways to similar and more strong results were presented, in particular, an approach 
based on general properties of t he L iouville evolution operator and "generalized 
fluctuation-dissipation relations" 15 16]. 
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